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A  Comprehensive  Description  of  the 
Mission  Sensor  Microwave  Imager  (SSM/I) 
Environmental  Parameter  Extraction  Algorithm 


I.  INTRODUCTION 

The  Mission  Sensor  Microwave  Imager  (SSM/I)  instrument  and  associated 
data  processing  software  are  being  acquired  as  part  of  the  joint  Navy/Air 
Force  project.  It  is  a  passive  microwave  radiometric  system  developed  by 
the  Hughes  Aircraft  Company  (HAC)  under  the  direction  of  the  Navy  Space 
Systems  Activity  (NSSA)  and  the  Air  Force  Space  Division  to  be  flown  on  the 
Defense  Meteorological  Satellite  Program  (DMSP)  operational  spacecraft  as  an 
near  all  weather  oceanographic  and  meteorological  sensor. 

The  SSM/I  is  a  seven-channel,  four-frequency,  linearly  polarized, 
passive  microwave  radiometric  system.  The  instrument  measures  atmospheric/ 
surface  brightness  temperatures  at  19.3,  22.2,  37.0  and  85.5  GHz  (1).  These 
data  will  be  processed  by  the  environmental  parameter  extraction  (EPE) 
algorithm  in  place  at  the  Fleet  Numerical  Oceanography  Center  (FN0C)  and  the 
Air  Force  Global  Weather  Center  (AFGWC)  to  obtain  near  real  time 
precipitation  maps,  sea  ice  morphology,  marine  surface  wind  speed,  columnar 
integrated  liquid  water,  and  soil  moisture  percentage.  These  products  will 
be  distributed  to  Navy  and  Air  Force  DMSP  operational  sites  to  satisfy  their 
unique  mission  requirements.  They  will  also  be  made  available  to  the 
general  scientific  and  industrial  communities  through  the  shared  METSAT  data 
archival  agreement  between  the  National  Oceanic  and  Atmospheric 
Administration  (NOAA)  and  the  Department  of  Defense  (DOD). 

The  SSM/I  data  processing  algorithm  was  developed  at  Environmental 
Research  A  Technology,  Inc.  (ERT)  under  a  sub-contract  from  HAC.  The 
algorithm  is  composed  of  four  modules  (2).  The  first  module,  SMISDP, 

Ingests  and  processes  raw  satellite  data  to  produce  earth-located  brightness 
temperatures  (Tg*s).  The  Tg»g  are  then  processed  through  the  second  module, 
Munucript  approved  August  10,  1983. 
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SMIEPE,  which  is  the  environmental  parameter  extraction  algorithm,  to  pro¬ 
duce  estimates  of  ocean,  land,  ice,  and  atmospheric  parameters.  The  third 
module,  SHIVER,  verifies  the  retrieved  parameters  against  pre-selected 
climatology  and  ground  truth  fields.  The  last  module,  SMIUPD,  updates  the 
coefficients  used  in  the  parameter  extraction  algorithm  when  the 
verification  results  meet  certain  pre-determined  criteria.  Large  volumes  of 
documentation  of  the  SSM/I  data  processing  algorithm  software  and  data  files 
exist  (e.g.,  3,4,5). 

The  environmental  parameter  extraction  module  retrieves  the 
environmental  parameters  using  the  brightness  temperatures  as  independent 
variables  in  linear  regression  equations.  The  regression  coefficients, 
which  reside  in  a  data  file  as  part  of  the  software  system,  are  determined 
using  geophysical  models,  radiative  transfer  models,  an  inversion  algorithm, 
and  climatology,  which  are  not  part  of  the  SSM/I  software.  The 
documentation  of  these  important  components  is  scanty  and  often  outdated 
(2, 6, 7,8) . 

The  purpose  of  this  report  is  to  present  a  comprehensive  digest  of  the 
EPE  algorithm  and  the  related  models  for  those  interested  in  the  SSM/I  and 
its  applications.  The  substance  for  this  report  has  been  taken  from  all  the 
previously  mentioned  references  and  a  software  listing  of  the  geophysical 
models  provided  by  Dr.  Kenneth  Hardy  (9),  the  former  leader  of  ERT's  SSM/I 
team.  This  report  Includes  an  overview  of  the  EPE  algorithm,  the  geo¬ 
physical  models  involved,  and  the  coefficients  and  important  criteria  used 
in  the  development  of  the  retrievals  of  environmental  parameters. 

II.  THE  DEVELOPMENT  OF  THE  SSM/I  EPE  ALGORITHM 

II. A  The  D  Matrix  Approach 

The  EPE  algorithm,  except  for  ice  morphology,  is  based  on  the  so-called 
D-matrlx  approach.  It  is  a  statistical  method  which  chooses  the  most 
probable  atmospheric  and  surface  properties  that  produce  the  set  of  measured 
brightness  temperatures  (Tg's).  The  formalism  for  the  statistical  approach 
starts  with  the  assumption  that  there  exists  some  linear  combination  of  T0's 
which  will  provide  information  about  the  geophysical  parameters  in  question. 


-  • 


That  is  to  say: 

n 

pi*  *  £  DIi  Ti»  where  i  ■  1,  - ,  m 

J*1  J  J 

or  in  vector  form 

P*  «  D*  T,  (II. 1) 

• 

where  P.  is  the  estimate  of  the  i-th  parameter  P . ,  T.  is  the  T_  of  the  j-th 

*  f  ^  J  ® 

radiometric  channel  and  D^j  is  the  (i,j)-th  element  of  the  D  matrix. 

The  matrix  used  in  the  SSM/I  EPE  algorithm  is  "tuned"  to  the  average 
value  of  each  parameter.  This  is  accomplished  by  transforming  the  vector  T 

into 

♦  (T)  ■  (1,  Tx  -  - ,  Tr  -  Tr).  (II. 2) 

The  effect  of  this  transformation  is  to  add  a  column  to  the  D  matrix  which 
contains  the  ensemble  averages  of  the  parameters. 

P*  «  D  ♦  (T),  (II. 3) 

where  £  is  the  new  D  matrix  with  dimensions  (m,  n+1). 

The  elements  of  the  D  matrix  are  defined  by  minimizing  the  mean  square 
error  between  the  predicted  and  the  actual  values  of  the  parameter  P^.  It 

can  he  shown  that 

■£  (P,  ♦)#£-1  (♦,  ♦),  (II.1*) 

where  £  (P,  ♦)  is  the  correlation  matrix  between  P  and  ♦  and  Cf1  (♦,  ♦)  is 

the  inverse  of  the  auto-correlation  matrix  of  ♦ . 

The  development  of  the  D  matrix  is  demonstrated  in  Figure  II. 1. 
Radiosonde  records  and  climatology  are  compiled  to  define  the  geophysical 


system  including  the  atmosphere  and  the  earth  surface  for  radiative  transfer 
calculations.  A  list  of  the  desired  environmental  parameters  is  given  in 
Table  XI. 1.  The  atmosphere  and  surface  data  are  used  to  simulate  Tg's  and 
thus  £(T)  through  a  geophysical-radiative  model.  .Half  of  the  P  and  d;  arrays 
are  used  to  oaleulate  correlation  matrices  and  the  D  matrix  (See  equation 

II. A).  The  other  half  of  the  $  array  is  used  to  test  the  validity  of  the  D 

# 

matrix  by  comparing  the  ground  truth  P  to  the  estimated  values,  P  which  are 
calculated  through  equation  II. 3.  The  D  matrices  installed  at  FNOC  and 
AFGIfC  have  been  tested  with  this  procedure. 

II. B  Separate  Schemes  for  Each  Climate  Zone 

The  assumption  of  linearity  between  P  and  T  is  not  satisfied  over  the 
entire  earth.  Instead  of  adding  higher  order  terms  in  the  regression 
equations,  a  number  of  climate  zones  have  been  selected  for  the  SSM/I  EPE 
algorithm.  For  each  of  these  zones,  linearity  is  assumed.  The  climate 
zones  are  defined  in  Table  II. 2.  AD  matrix  is  developed  for  each  of  the 
climate  zones.  The  D  matrices  of  the  transition  zones,  are  averages  of 
those  of  neighboring  zones. 

Due  to  the  limitation  of  computing  facilities  at  AFGWC,  the  maximum 
number  of  channels  used  for  the  retrieval  of  each  parameter  is  limited  to 
four.  The  channel  selections  for  all  the  parameters  are  shown  in  Table 
II. 1.  For  example,  the  regression  equation  for  the  sea  surface  wind  speed 
la 

SW  *  do  ♦  d1  Tb(19H>  ♦  d2  Tg(22V)  +  d3  Tg(37V)  ♦  d^  Tg(37H), 
where  d  .  . ...,d„  are  D  matrix  coefficients.  The  value  of  the  coefficients 

O  *1 

for  each  environmental  parameter  for  all  climate  zones  will  be  described  in 
Section  IV. 

II. C  Piecewise  Scheme 

Even  though  all  the  environmental  parameters  for  each  climate  zone  car 
be  simultaneously  extracted,  the  results  may  not  e  meaningful.  For 
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TABLE  II. 2  Designation  of  SSM/I  Climate  Zones 
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Ocean 
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13 
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Ocean 
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14 
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Sept 

Oct 
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5 

5 

9 

8 

6 

6 

6 

10 

11 

11 

11 

13 

13 

12 

12 

12 

14 

14 

1:  tropical  warm/land 
2:  tropical  warm/ocean 
3:  tropical  cool/land 
4:  tropical  cool/ocean 
5:  mid  lat  spring-fall/land 
6:  mid  lat  spring-fall/ocean 
7:  mid  lat  summer/land 
8:  mid  lat  summer/ocean 
9:  mid  lat  winter/land 
10:  mid  lat  winter/ocean 
11:  arctic  cool/land 
12:  arctic  cool/ocean 
13:  arctic  cold/land 
14:  arctic  cold/ocean 
15:  transition  zones 


A.  Lower .Latitude  transition  zones  (LLTS)  are  between  tropics  and  mid  lat. 

a.  LLTS  warm/land 

b.  LLTS  warm/ ocean 

c.  LLTS  cool/land 

d.  LLTS  cool/ocean 


B.  Upper  Latutide  transition  zones  (ULTS)  are  between  mid  lat.  and  arctic 


a.  ULTS  cool/land 

b.  ULTS  cool/ocean 

c.  ULTS  cold/land 

d.  ULTS  cold/ocean 
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example,  if  there  is  heavy  rain  over  land,  the  brightness  temperatures  will 
not  be  sensitive  to  the  cloud  water  and  soil  moisture  conditions. 

Similarly,  when  there  is  heavy  rain  over  ocean,  cloud  water  and  sea  surface 
wind  retrievals  would  not  be  reliable.  The  SSM/I  EPE  algorithm  employs  the 
so-called  piecewise  scheme,  to  determine  which  parameters  to  retrieve  based 
on  the  values  of  certain  Tg  observations.  The  piecewise  schemes  for  ocean 
and  land  are  depicted  in  Figures  II. 2  and  II. 3.  The  choice  of  parameters  to 
retrieve  is  mainly  based  on  rainfall.  The  criteria  for  determining  whether 
there  is  no  rain,  light  rain,  or  heavy  rain  are  different  for  each  climate 
zone.  They  are  listed  in  Section  IV. 

The  retrieval  of  ice  parameters,  including  concentration  and  age,  is 
accomplished  through  a  deterministic  approach.  The  ice  algorithm  is 
developed  directly  from  physical  relationships  rather  than  through 
regression.  The  description  of  the  physics  of  ice  morphology  and  retrieval 
will  be  presented  in  Sections  III  and  VI. 

III.  THE  GEOPHYSICAL-RADIATIVE  MODEL 

The  radiative  transfer  theory  is  used  to  simulate  the  brightness 
temperature  obtained  by  a  remote  sensing  platform,  such  as  a  satellite, 
observing  the  earth-atmosphere  system.  A  geophysical  model  which  specifies 
the  physical  and  electromagnetic  properties  of  the  atmospheric  and  the 
earth’s  surface  provides  the  necessary  input  parameters  to  the  radiative 
transfer  theory. 

III. A  The  Radiative  Transfer  Model 

For  microwave  radiation,  the  important  atmospheric  parameters  are  water 
vapor,  oxygen  and  liquid  water,  both  in  the  form  of  cloud  water  and  rain. 

The  form  of  the  liquid  water  governs  the  method  of  calculation  for  the 
extinction  coefficient  and  the  solution  to  the  radiative  transfer  equation. 
When  liquid  water  exists  in  small  droplets  such  as  found  in  clouds,  the 
attenuation  is  predominantly  by  absorption  and  calculations  are  relatively 
straightforward .  But  in  the  case  of  rain,  water  drops  are  so  large  that 
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Figure  II. 3  Summary  of  piecewise  algorithm  over  land 


scattering  effects  must  be  considered.  The  Mie  theory  is  employed  to 
calculate  the  extinction  coefficient  including  both  absorption  and 
scattering  effects.  An  interative  algorithm  is  needed  to  solve  the 
radiative  transfer  equation. 

The  radiative  transfer  equation,  with  the  Rayleigh-Jeans  long 
wavelength  approximation,  for  the  intensity  of  thermal  radiation  from  a 
blackbody  can  be  written  as 


dTg(Z,0) 


Tb(Z,0)  -  T(Z) , 


(III. 1) 


where  Tg  is  the  brightness  temperature,  a  function  of  height,  Z,  and  of  look 
angle  9.  T  is  physical  temperature  and  r  is  the  opacity.  In  the  absence  of 
scattering,  the  Tg  sensed  by  a  satellite  can  be  expressed  as  (9) 


Tg(Z.g)  =  Tg  (Z,0)  +  [(1-R)Tg  +  R  Tb  ]  e  _T  , 


(III. 2) 


where  Z  is  the  height  of  the  satellite,  R  is  the  effective  surface 
reflectivity,  r  is  the  total  opacity  of  the  atmosphere  along  the  line  of 
sight  (nepers)  and  T_  is  the  surface  temperature  (K).  The  quantities  Tn 

1 

and  Tn  are  proportional  to  the  upward  emission  from  the  atmosphere  and  the 
B2 

downward  emission  from  the  atmosphere  plus  attenuated  sky  background 
radiation  and  are  given  by 
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(III. 3) 


T-  *  T 


e-T  ♦  f 


(-fy  (z1)  seefl  dz') 


T(z)  (z)  e 


seed  dz. 


(III. 4) 


where 


1 


y(z)  sec$  dz  and  “Kz)  is  the  total  opacity  at 


height  z,  representing  the  sum  of  contributions  from  water  vapor,  oxygen, 
and  cloud  water  droplets. 

Equations  III. 2  through  III. 4  are  valid  for  all  non-precipitating 
atmospheres.  When  precipitation  is  present  multiple  scattering  effects  are 
coupled  with  absorption.  The  Mie  theory  is  based  on  the  diffraction  of  a 
plane  monochromatic  wave  by  a  sphere  with  a  homogeneous,  complex  index  of 
refraction.  By  examining  the  dissipation  in  the  sphere  and  the  scattered 
wave  produced,  it  provides  the  information  necessary  for  the  calculation  of 
the  extinction  coefficient  and  single  scattering  albedo. 

The  Mie  efficiency  factors  for  scattering,  absorption  and  extinction 
are  defined  as  the  ratio  of  the  actual  to  the  geometric  cross-sections. 


Qj  •  V" 


J  =  S  (scattering) 

A  (absorption) 

E  (extinction) , 


(III. 5) 


where  X  is  the  actual  cross-section  and  r  is  the  raindrop  radius.  By 
definition  QE  s  QA  +  Qg.  The  single  scattering  albedo  is 


•o  =  VQA 


(III. 6) 


The  extinction  coefficient  is  determined  by  integrating  the  efficiency 
factor  over  the  drop  size  distribution  encountered  in  the  rain. 


=  N(r)  Q£  (n  (A.T),  A,r)  r  dr, 


(III. 7) 


where  N(r)  is  the  number  of  drops  of  radius  r  per  unit  volume,  A  is  the 
observational  wavelength,  rf  is  the  index  of  retraction  and  T  is  the 
temperature.  The  radiative  transfer  equation  for  the  scattering  medium  can 


be  written  as 


dT-(  T\M) 

m - - -  =  T  (r',M)  -  T_(t\m)  , 

d  »  B  S 


(III. 8) 
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where  a  is  the  cosine  of  the  incidence  angle,  r*  =  r/fj  and  Tg  is  the 


source  function 


Ts  *  n  --o 


wo(T,)  f 1 

(  r* )  ]  T  (t',M) - -  Jy  p(M.M')  Tg  (t\m»)  dM\ 


(III. 9) 


where  P  is  the  phase  function.  The  integral  equation  for  T_  can  be  written 


as 


-  T  U>  (  T  )  /*T 
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r  /*r  M  i 

r«.r)  (1-R)  Tg  ♦  KJ  Ts(r)  e  jjL  .(III. 


10) 


where  *  is  the  total  optical  depth  of  the  atmosphere  and  E^  is 


v  / w \  1  ("*x/m  )  n-2  . 

E  (x)  a  f  e  M  <V. 

n  o 


(III. 11) 


A  variational-iterative  approach  (10,11)  is  used  for  the  solution  of  the 
source  function  Tg  defined  in  equation  (III. 10). 


III.B  The  Atmosphere  Model 


The  physical  electromagnetic  properties  of  water  vapor,  oxygen,  and 
liquid  water  as  defined  in  the  ERT  model  are  described  in  this  sub-section. 


III.B. 1  Water  Vapor 

The  extinction  coefficient  due  to  water  vapor  as  described  by  Gaut  (12) 
is  used.  Laboratory  measurements  and  theoretically  generated  data  were  used 
to  examine  and  modify  the  expression  formulated  by  Van  Vleck  and  Weisskopf 
(13).  Only  the  absorption  line  centered  at  22.235  GHz  is  used  to  calculate 
the  extinction  coefficients  for  all  but  the  highest  frequency,  i.e.,  85.5 
GHz,  of  the  SSM/I  instrument.  The  extinction  due  to  water  vapor  is  composed 
of  two  distinct  portions;  the  resonance  and  non-resonance  components.  Thi 


former  is  the  rapidly  changihg  contribution  from  the  nearest  water  vapor 
absorption  line  while  the  latter  is  due  to  the  slowly  varying  absorption 
wings  of  more  distant  lines.)  \ 


*  tWfBES  (2)  ♦  YWfN0N(^) 


.jk,  pM<Z)  V2  T(z)-2-5le  WWW) 


(v-v  )  +  Av 

o 


2  2 
(v+v  )  +  Av 

o 


+  5.k2  pw(z)v2  T(z)"1,5av  !  X  100  , 


(III. 12) 


where  yw(z)  is  the  absorption  due  to  water  vapor  at  height  z  (nepers/m), 
v  is  frequency  (GHz),  vq  is  the  reference  frequency  (22.235  GHz)  and  k1t 
are  constants. 

k1  =  3.5175  x  10~3;  k2  =  5.0920  x  10-9. 

The  half-width  of  the  absorption  line,  Av,  is  defined  as 

Av  *  2.62  (P(z)/PQ)  (T(z)/To)"°*625  [1  +  0.15  pw(z>  T(z)/P(z)],  (III. 13) 

where  P(z)  is  pressure  at  height  z  (mb),  P  is  1013.25(mb),  and  T  is  318.0 

o  o 

(K). 

For  85.5  GHz,  eight  water  vapor  rotational  spectral  lines  are  included 
in  the  calculation  of  the  absorption  coefficient.  The  reference  frequencies 
of  these  lines  and  other  significant  parameters  are  listed  below. 


i 

V 

0 

Stn(i) 

C^i) 

c2(i> 

t  (i) 
exp 

t1(i> 

t2(i) 

Sp(i) 

1 

22.235 

.0549 

2.7057 

.01976 

.626 

446.39 

447. 17 

3 

2 

183.310 

.1015 

2.8800 

.01906 

.649 

136.15 

142.30 

1 

3 

323. 159 

.0870 

2.2956 

.01952 

.420 

1283.02 

1293.80 

3 

4 

323.758 

.0891 

2.7876 

.02050 

.619 

315.70 

326.50 

1 

5 

377.418 

.1224 

2.8440 

.02102 

.630 

212. 12 

224.71 

3 

6 

389.709 

.0680 

2.1060 

.02035 

.330 

1525.31 

1538.31 

1 

7 

435.874 

.0820 

1.5000 

.01976 

.290 

1045.  14 

1059.68 

1 

8 

437.673 

.0987 

1.7700 

.02253 

.360 

742.18 

756.18 

3 

The  resonance  portion  of  the  absorption  coefficient,  ^(z) ,  is  given 
y  (z>  *  £  *>5.5  «,(z>  Sp(i)  Stn(i)|"- - -1 

w.Res  1.1  *  Lt(z)1-5  J 

(1>1  [  T7T7?  . 


t'2  -  -2< 

o 


where  Sp(i)  is  the  spin  mode  and  Stn(i)  is  the  line  strength. 


A* 


■e,<u?r)W 


W‘> 


El  +  C2(i)  ^(z)  T(z)/P(z)] 


where  t  is  the  exponent  for  the  temperature  term  and 
exp 


tm  =  exp 


r  t(i)-|  r  t<i)i 

[ -1' •' 43897  -JUT}-  exp  [-1.43879-^J 


[t,(i>  -  t  (i)] 

if  -1.43879  — - -  >0.1, 


T(z) 


tm(i)  •  ‘'o<1)  ^TTzT  exp 


[-U3879(?fir)] 


Ct  ( i)  -  t  (i) ] 

if  -1.43879  T(z)  g - <0.1. 

The  non-resonance  term  is  defined  as 


,-9 


VNon(l>  *  ’-2375  *  10  « 


<-m) 


1.5  + 


8 

l 

i  =  1 


Vt  (i) 
Z-r  exp 


and  the  total  absorption  due  to  water  vapor  is  then 


by 


(III. 14) 


(III. 15) 


(III. 16a) 


(III. 16b) 


(III. 17) 


r  (z)  *  r  „-_(z)  ♦  r.  ..._(z) 


(III. 18) 


III.B.2  Oxygen 


The  Van  Vleck  and  Weisskopf  (13)  formulation  for  the  absorption  line 
shape  is  adopted  for  the  calculation  of  absorptidn  due  to  oxygen. 


-3 


yQ(z)  =  C1  P(z)  (760/1013.25)  T(z)  -»2 

1  s  2  |fT[ i(2i+3)/( i+1 ) ]+f”[ ( i+1 ) (2i-1 )/i]+2f  [ ( i2+i+1 ) (2i+1 )/( i( i+1 ) )  1 
1,45,2 11  1  °  7 

e-2. 06844  i( i+1 )/T(z)  , 

where  f  *  =  A*  {  1/[(0+(i)  -  v)2  +  Ak2J  +  1./[(0+(i)  +  „  )2+  i 


(III. 19) 


and 


f.  =  A*  1  / C (0~( i)  -•>)*■  +  Ai/  ']  +  1./[0"(i)  +  v  )2  +A^2]  , 


(III. 20) 


«  P(z)(0.21  +  0.78  {) 


0.85 


(III. 21) 


where  s  0.25  if  P(z)  >  365.51  mb 
=  0.75  P(z)  <  25.5  mb 

=  0.25  ♦  0.5  £n^^  if  25.5  <  P(z)  <  365.51; 

f  =  Av/(i'2  +  Av2); 

®  .3 

C1  =  0.61576  x  10  J  J 

=  (1.95  x  760/1013.25)  x  10“3, 

and  0+(i)  and  0  (i)  are  line  shape  parameters  for  the  absorption  lines. 
They  are  defined  by  Rosenkrantz  (15)  and  are  listed  in  Table  III. 1 . 
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The  resonance  portion  of  the  absorption  coefficient,  y  (z),  is 

w 


8  I*  2  1 

y  <z>  =  £  45.5  .  (z)  SpU)  Stn(l)  - _ I 

1=1  Li(z)’-5  J 

(  v 2  -  v 2)2  +  Hv2A v2  J  m 


[r2  -  v2(i)3 
o 


where  Sp(i)  is  the  spin  mode  and  Stn(i)  is  the  line  strength. 


t  (i) 

A(zA/t  \exp 

A.s  El  +  C2(i)  pw(z)  T(z)/P(z)]  , 

where  t  is  the  exponent  for  the  temperature  term  and 

6*  p 

r  t(i)-j  r  tun 

tm  =  exp  -1.43897  -  exp  1-1 . 43879  ~fuyj 

tt  (i)  -  t_(i)] 

if  -1 .43879  — - - yuT - 2  0.1, 


ct  (i>  - 1  (i)  ] 

if  -1.43879  — '  T(z)  - - <0.1. 

The  non-resonance  term  is  defined  as 


V*on(l)  =  '-2375  x  ,0‘9 


1.5  ♦ 


8 

l 

i  =  1 


EW‘> 


and  the  total  absorption  due  to  water  vapor  is  then 


given  by 


(III. 14) 


(III. 15) 


(III. 16a) 


(III. 16b) 


m2,  (III. 17) 


*«<Z)  s  P.«(Z>  +  Won(z)- 


(III. 18) 


III. ft. 2  Oxygen 


The  Van  Vleck  and  Irfeisskopf  (13)  formulation  for  the  absorption  line 
shape  is  adopted  for  the  calculation  of  absorptidn  due  to  oxygen. 

~3  2 

y  (z)  s  C,  P(z)  (760/1013.25)  T(z)  ^ 
o  1 

is  53  i  (2i+3)/(  i+1 )  3+f7[  ( i+1 )  (2i 

1 ,45*2  "  1 

-2.06844  i( i+1  )/T(z)  , 

e 

where  f*  =  to  {  1/[  (0+(i)  -  v  )2  +Aw2]  +  1./[(0+(i)  +  „  )2+  A^2]J 

and 

f”  =  Ak  1/[(0~(i)  -»)2  +  to2}  +  1./[0~(i)  +  v)2  +  to21  .  (III. 20) 

fiv  =  «-P(z) (0.21  +  0.78  ,  (III. 21) 

where  =  0.25  if  P(z)  >  365.51  mb 
=  0.75  P(z)  <  25.5  mb 

=  0.25  +  0.5  if  25.5  <  P(z)  <  365.51; 

f  -=  to/(2  *  to2)', 

°  -3 

C1  =  0.61576  x  10  3  ; 

=  (1.95  x  760/1013.25)  x  10-3, 

and  0+(i)  and  0-(i)  are  line  shape  parameters  for  the  absorption  lines. 

They  are  defined  by  Rosenkrantz  (15)  and  are  listed  in  Table  XI 1.1. 


-1)/i]+2fo[(i2+i+1)(2i+1)/(i(i+1))|. 

(III. 19) 
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v-V*- 


TABLE  III. 1  Lin®  Shape  Parameters  for  the  Absorption  Due  to  Oxygen 


if 


i 

0*(i) 

0~(  i) 

1 

56.2648 

118.7505 

3 

58.4466 

62.4863 

5 

59.5910 

60.4863 

7 

60.4348 

59.1642 

9 

61.1506 

58.3239 

11 

61.8002 

57.6125 

13 

62.4112 

56.9682 

15 

62.9980 

56.3634 

17 

63.5685 

55.7839 

19 

64.1272 

55.2214 

21 

64.6779 

54.6728 

23 

65.2240 

54.1294 

25 

65.7626 

53.5960 

27 

66.2978 

53.0695 

29 

66.8313 

52.5458 

31 

67.3627 

52.0259 

33 

67.8923 

51.5091 

35 

68.4205 

50.9949 

37 

68.9478 

50.4830 

39 

69.4741 

49.9730 

41 

70.0000 

49.4648 

43 

70.5249 

48.9582 

45 

71.0497 

48.4530 

XII. B. 3  Liquid  Water 


a.  Cloud  droplets 

The  absorption  coefficients  of  cloud  droplets  for  all  frequencies 
exoept  85.5  GHz  are  calculated  based  on  the  formulation  by  Staelin,  et.  al . 
(14).  When  there  is  no  precipitation. 


z)  s  Cw(z)*10 


[0.0122*(291-T(z))  -4.0]/(30A)2 


(III. 22) 


where  Cw(z)  is  the  cloud  water  density  at  height  z. 

For  the  85.5  GHz,  the  Raleigh-Jeans  approximation  is  used  for  the 
calculation  of  the  absorption  coefficient 


yr(z)  =  1.88685  x  10“3  C  (z)*I-/(30A> ,  (III. 23) 

U  W  I 

2  2 

where  If  is  the  imaginary  part  of  (-m  +1)/(m  +2)  and  m  is  the  complex  index 
of  the  refraction  of  water. 

The  absorption  of  ice  cloud  particles  is  calculated  according  to  the 
following  formula,  assuming  the  temperature  is  below  freezing,  or 
T(z)<273. 16  K. 


T,(z)  =  [0.2302585-C  (z)*10 
v.  w 


-(5.11  +  0.145V273.16-T  ) 


]/(30/*). 


(III. 24) 


b.  Precipitation 

Precipitation  Is  the  primary  contributor  to  atmospheric  attenuation 
at  microwave  wavelengths.  The  attenuation  results  from  both  absorption  and 
scattering  by  the  hydrometeors.  The  magnitude  of  these  processes  depends 
upon  wavelength,  drop  size  distribution,  and  precipitation  layer  thickness. 

In  the  case  of  light  rain,  when  it  is  permissible  to  neglect  the 
effect  of  multiple  scattering,  the  attenuation  due  to  precipitation  can  be 
calculated  according  to  equation  III. 12  and  III. 18.  The  liquid  water  con- 
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-1  j 


tent  and  the  rain  rate  are  assumed  to  be  related  by  the  Marshall-Palmer 
distribution  function 


1.28  x  10“10  x  • 


A$0  RR*21\ 
\  -3.67  / 


(111.25) 


where  Lw  is  liquid  water  content  (gm/in  )  and  RR  is  rain  rate  (mm/hr.) 

For  heavier  rain,  the  multiple  scattering  effect  will  be  included 
through  the  algorithm  described  in  sub-section  III. A,  equations  III. 5  thru 
III. 11.  One  of  the  essential  elements  for  the  algorithm  is  the  drop  size 
distribution.  The  empirically  observed  spectra  of  Laws  and  Parson  (l6)  is 
used  to  fit  the  DeirmendJ ian  drop  size  distribution  for  the  SSM/I  rain 
algorithm.  The  DeirmendJ ian  distribution  is  defined  as 


n(r)  ®  Ar^  e  ^~Br  ^  ,  (111.26) 

where  n(r)  is  the  drop  size  distribution  (cm~3Uni  ),  r  is  the  radius  of 
raindrop  (tun),  and  A  and  B  are  scale  parameters. 


[(^  +  h)/c2\ 


where  M  is  the  total  liquid  water  content  (gm/ni  ), 


(kt  3 

Mw=  pwj  — *  r  n(r)dr. 


(III. 27) 


(III. 28) 


@-o2 
ro 


(111.29) 


and  rQ  is  the  mode  radius  (Mm).  The  shape  parameters  are  and  C2.  The 
former  of  these  affects  the  distribution  of  the  smaller  radii  while  the 
latter  affects  the  larger  radii. 
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The  analytical  results  from  fitting  the  Laws-Parson  empirical  spectrum 

are: 

M  =  0.636  (RR'®®1)  (gm/cm^),  where  RR  is  rain  rate  (mm/hr), 

r  =  225  +  9.16  (RR*881),  (r  )  =  375  („m), 

o  o  max 

C1  =  4.0,  (III. 30) 

C_  =  0.70  -  0.00458  (RR*881),  (CD  .  =  .625. 

k  i.  mxn 

This  formulation  relates  the  rain  rate  to  mode  radius  and  liquid  water 
density.  It  parameterizes  the  distribution  as  a  function  of  rain  rate  only 
and  thus  simplifies  the  solution  to  the  Deirmendjian  distribution.  On  the 
other  hand,  it  still  retains  the  sensitivity  of  scattering  to  drop  size, 
which  would  be  lost  if  the  Marshall-Palmer  distribution  is  used.  Figure 
III.1  shows  the  good  agreement  between  the  Deirmendjian  model  and  the  Laws 
and  Parson  distributions. 

For  weak  storms  with  low  rain  rate  (RR  <  8  mm/hr),  the  ceiling  for  the 
rain  layer  is  assumed  to  be  the  0  C  isotherm.  For  more  intensive  convective 
storms,  the  ceiling  is  usually  higher  with  a  super-cooled  water  layer  on 
top.  The  thickness  of  the  layer,  Zj  ,  depends  on  the  intensity  of 
convection.  In  the  ERT  algorithm,  the  thickness  of  the  supercooled  layer, 
AZ,  for  rain  rates  in  excess  of  8  mm/hr,  is  given  as 

AZ  =  (RR  -  8.0)  Z  :  AZ  <AZm.  (III. 3D 

inc  max 

where  Z.  and  AZ  are  determined  from  climatology  (17).  The  values  used 
inc  max 

are  given  below: 
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I  2  3 

DROP  RADIUS  (mm) 


The  matching  of  the  Deirmendjian  raindrop  size  distribution 
used  in  model  computations  to  the  observed  Laws  and  Parsons 
distributions. 


Z.  (m/hr) 

me 

mm 

AZ  max  (m) 

Tropics  Land 

50 

800 

Ocean 

50 

800 

Mid  Lat  Land 

100 

1600 

Ocean 

50 

800 

Arctic  Land 

25 

400 

Ocean 

25 

400 

III.C  The  Surface  Model 

III.C.1  Ocean  Surface 

The  calculation  of  the  emissivity  of  a  smooth  water  surface  is 
relatively  straightforward.  The  dielectric  properties  of  sea  water,  derived 
from  the  measurements  of  Lane  and  Saxton  (18)  and  expressed  in  analytical 
form  by  Chang  and  Wilheit  (19)  are  used  for  the  SSM/I  algorithm.  The 
Fresnel  equations  for  a  plane  dielectric  interface  are  used  to  calculate  the 
emissivity  of  a  smooth  surface  for  a  given  view  angle  and  polarization. 

Wind  driven  waves  and  foam  on  the  ocean  surface  both  significantly 
alter  the  microwave  emissivity  of  the  ocean  surface.  Wind  speed  is  highly 
variable  both  horizontally  and  vertically.  Marine  wind  speed  measured  from 
a  ship  is  usually  referenced  to  a  standard  height  of  20  meters.  However, 
the  wind  speed  which  directly  relates  to  roughness  and  foam  is  the  friction 
velocity  at  the  ocean-air  interface.  A  relationship  developed  by  Cardone 
(20)  is  used  to  extrapolate  from  the  standard  20  meter  height. 

The  ocean  surface  model  developed  by  Wilheit  (21)  is  adopted  in  the 
33M/I  algorithm.  The  roughness  effect  is  modeled  after  Cox  and  Munk  (22) 
tremtiag  the  sea  surface  as  a  collection  of  plane  fp'vts  large  compared  to 

the  observational  wavelength.  The  variance  in  se.  .  'ace  slope  is  defined 

» 

am  a  function  of  win'*  speed  and  then  used  Ifor  calculation  of  emissivity 
assuming  a  Gaussian  lope  distribution  fon  the  face  i  using  the  Fresnel 

relations.  The  deri'  i  wind  speed  dependence  of  the  s» rface  slope  variance, 

2 

•  ,  is  given  by 


#2(f)  =  (0.3  +  0, 02f )  (0.003  ♦  0. 48w) 


for  f  <  35  GHz; 
for  f  >  35  GHz, 


(III. 32) 


B 


•>.i 


3 

$ 


1 

1 

i 
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=  0.003  +  0.48w 

where  f  is  frequency  in  gigahertz  (GHz)  and  w  is ‘the  wind  speed  (m/sec)  at 
20  meter  height. 

Foam  is  treated  as  partially  obscuring  the  surface  in  a  manner  indepen¬ 
dent  of  polarization  but  dependent  upon  frequency.  The  foam  fraction,  K,  is 
defined  as 

-f/7  5 

K  =  0.006  (w  -  7)  O  -  e  ' for  w  >  7  m/sec 

=0  w  <  7  m/sec.  (III. 33) 

This  equation  is  based  on  Nordberg’s  (23)  observation  that  the 
brightness  temperature  increases  linearly  with  wind  speed  exceeding  7  m/sec, 
and  that  no  foam  forms  below  that  wind  speed. 

III.C.2  Land  Surface 

The  description  of  the  land  surface  is  extremely  complex  due  to  the 
many  types  of  surfaces  and  the  variation  of  the  physical  characteristics 
within  each  surface  type.  At  present,  the  only  land  surface  parameter 
retrieved  by  the  SSM/I  algorithm  is  the  soil  moisture  for  open  land  areas. 
Criteria  have  been  developed  to  separate  other  land  surface  types  such  as 
forest  and  lakes  from  the  open  land  areas. 

The  intensity  of  radiation  from  3oil  depends  on  the  local  dielectric 
constant  and  the  physical  temperature  of  the  soil.  Moisture  produces  a 
marked  increase  in  both  the  real  and  the  imaginary  parts  of  the  dielectric 
constant  of  soil,  leading  to  a  decrease  in  the  soil  emissivity.  Experi¬ 
mental  observations  and  theoretical  calculations  indicate  that  the 
emissivity  of  soil  at  microwave  frequencies  can  range  from  >0.9  for  dry 
soils  to  £0.6  for  very  moist  soils. 

A  generalized  incoherent  layered  surface  model  (24),  was  used  to 
simulate  expected  SSM/I  soil  moisture  brightness  temperature  signatures. 
Assume  that  there  are  N  soil  layers,  with  the  atmospheric  layer  above  and  an 
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infinitely  deep  layer,  with  constant  temperature  and  moisture  content, 
below.  The  brightness  temperature  emerging  at  the  si  il  surface  for  the 
horizontal  polarization,  TgH,  is  given  by 


Tbh  =  p  T(i)  W(i)  T^CiJIl  -  e"s(l)  dz(i)]  [1  +  RH(i)  •  g(i)  dz(i)|  jni>- 

where  T(i)  is  the  temperature  of  the  i-th  layer  and  dz(i)  is  the  depth  of 
the  i-th  layer.  The  weighting  function  w(i),  is  defined  as 


w(i)  =  e 


i=l 

-  £g(m)  dz(m) 

m=l 


(III. 35) 


where  g(m)  =  k*v/c  •  A(m).  Here  c  is  the  speed  of  light 

and  A(m)  »  «j(m)/[2  B(m)l,  where  *j(m)  is  the  imaginary  part  of  the 

dielectric  constant  of  the  m-th  layer. 

/  2  1/2  >1/2 
B(m)  *  \  0.5  Y(m)  •  [l  +  (1  +  («I(m)/Y(m )T)  ]}  ,  (III. 36) 

2 

and  Y(m)  *  *D(m)  -  sin  •,  where  «  (m)  is  the  real  part  of  the  dielectric 
n  h 

constant  of  the  m-th  layer  and  6  is  the  incidence  angle.  The  term,  T^i), 
is  the  transmission  factor  of  the  i-th  layer  for  horizontal 
polarization  and  is  given  by 


T^U)  =  [1  -  Rgd)] . [1  -  yi-l)] 


(III. 37) 


where  Rg(i)  is  the  horizontal  extinction  coefficient  for  the  i-th  layer. 


V1’ 


C -  Cw(i  *  11 
Cw<1)  +  C«v(1  +  11 


,  where  C^d)  =  (B(i),  A(i)). 


(111.38) 


The  brightness  temperature  for  vertical  polarization,  Tgy.  emerging 
at  the  surface  is 

T  *  E  T(i)  W(i)  T  (i)  (1  -  e"g(l)  dz(l)]  [l  +  Ry(i)  g(i)  dz(i)].  (III. 39) 


Here  TRy(i)  is  the  transmission  factor  of  the  i-th  layer  for  vertical 

polarization; 

TRV(i)  “  (1  -  V1*  (1  "  V2)) . (1  '  Vi-l))’  (111.40) 

where  Ry(i)  is  the  vertical  extinction  coefficient  for  the  i-th  layer. 
e(i)Cw(i+l)  -  e (i+l)  Cwv(i) 

Vi)  "  - 

e(i)Cwy(i+l)  +  e (i+l)  C^i) 

The  dielectric  constant  is  a  function  of  frequency  and  soil  moisture. 

For  the  19*35  and  22.235  GHz  channels,  results  from  the  study  by  Wang  et.al. 

(25)  are  adopted.  Let  SM  be  soil  moisture  in  percentage.  The  real  part  of 
the  dielectric  constant,  eR,  is  given  by 

SM  <_  7  eR  =  2.6?  +  0.1  SM, 

7  <  SM  <_  19.25  eR  *  3.37  +  0.78  (SM  -  7),  (Ill. 1+2) 

19.25  <  SM  er  *  12.92  +  0.23  (SM  -  19-25), 

and  for  the  imagniary  part,  e^,  by 

SM  <  10  e  =  0.08  SM, 

10  <  SM  _<  19  Cj  =  0.80  +  0.57  (SM  -  10),  (111.43) 

19  <  SM  e  =  5.93  +  0.28  (SM  -  19). 

The  results  from  the  study  by  Geiger  and  Williams  (26)  are  adopted  for 

the  dielectric  constant  of  soil  at  37*0  and  85.5  GHz.  Here 

eR  »  (1/0.885) (-2.8426  +  1.0534  SM  -  0.033808  SM2  +  3.434  x  lO"*4  SM3) 

and  (III. 44) 

*  (1/0.885) (-1.2981  +  0.2925  SM  +  1.79  x  10“3SM2  -  1.465  x  10-USM3). 


2 

,  where  e(i)  =  (eR(i),  e^i)).  ( III  .4 1 ) 


At  the  higher  frequencies,  there  is  little  sensitivity  to  soil  moisture 
much  below  the  earth  surface.  It  suffices  to  calculate  the  reflection 
coefficient  through  the  Fresnel  coefficients  using  the  dielectric  constants 

specified  above. 

Choudhury  and  Schmugge  (27)  have  developed  a  simple  correction  factor 
to  the  reflectivity  to  take  the  surface  roughness  into  account.  The 
corrected  reflection  coefficient  is  given  by 

R'(9)  =  R(e)  e(‘h  cos  0),  (in. 46) 

where 

h  *  4o2(2*/x)2. 

2 

a  is  the  variance  of  the  surface  roughness,  X  is  the  observational 
wavelength  and  h  at  19.0  GHz  is  assumed  to  vary  between  0  and  0.6. 

III.C.3  Sea  Ice 


The  upwelling  brightness  temperature  of  a  scene  containing  sea  water 
and  various  amounts  of  sea  ice  is  a  function  of  the  ice  concentration,  ice 
emissivity,  the  physical  temperatures  of  the  ice  components,  and  the  amount 
of  water  vapor  and  liquid  water  in  the  atmosphere  (6,7).  The  ice  can  be 
composed  of  first  year  (FY),  multi-year  (MY)  or  first-year  thin  ice  (FT). 
The  vertically  polarized  up-welling  brightness  temperature  sensed  by  the 
S8M/I  radiometer  is  expressed  as 


To,,  *  e"T  [C*F(c  -  T_  +  (l-e  )T  )  +  C*M(e  T  +  (l-£  )T  ) 
BV  fv  f  fv  s  mv  m  mv  s 


*  C(l— f-M)UtvTt  *(l-«tv)Ta)  *  (l-CHc^T,  ♦  U-„)T,) 


(III. U6) 


where  r  =  total  atmospheric  opacity, 


C  s  fraction  of  sea  ice  (includes  all  types  within  field  of  view), 

F  =  fraction  of  FY  ice  relative  to  total  ice  present, 

M  =  fraction  of  MY  ice  relative  to  total  ice  present, 

«,  , «  , =  vertically  polarized  surface  emissivities  of 

Z  V  inv  WV 

FY,  MY  and  FT  ice, 

w  =  vertically  polarized  surface  emissivity  of  sea  water, 

Tf,T|n,Tt  =  surface  temperature  of  FY,  MY  and  FT  ice, 

Tw  s  surface  temperature  of  sea  water, 

T  s  incident  sky  temperature  at  the  surface  due  to  atmospheric 

3 

downward  self-emission,  and 

T  =  contribution  from  atmospheric  upward  self-smission . 

Sum 

An  approximate  expression  for  r  at  a  50°  earth  incidence  angle  for 
( *  a  6.6  to  37  GHz)  is 

r  =  L,  (III. 47) 

1C 

where  A  is  the  observational  wavelength  (cm)  and  L  is  the  integrated 
vertical  column  of  liquid  water  and  water  vapor  (cm).  An  expression  similar 
to  equation  III. 46  exists  for  horizontal  polarization.  Ideally  all  SSM/I 
frequencies  with  different  sea  ice  dependencies  can  be  used  to  solve  for 
C,  F,  and  M  simultaneously.  However,  limitations  at  the  AFGWC  computing 
facility  prevent  the  use  of  all  the  SSM/I  channels  in  a  single  algorithm.  A 
simpler  but  reliable  algorithm  was  developed. 

Results  from  a  study  on  ice  emissivities  (28)  are  presented  in  Table 
III. 2.  The  difference  between  the  vertical  and  horizontal  emissivities  at 
37  GHz  exhibits  significantly  less  dependence  on  ice  type  than  does  the 


27 


Emissivity  of  ice 


rr  < 

x  uj 

uj  co 

S5 

s< 


to 

X 

CO 

o 

o 

o 

o 

> 

CD 

o‘ 

o 

CD 

X 

00 

o 

o 

CD 

CD 

> 

0) 

o> 

o 

o 

ID 

X 

00 

C5 

o 

d 

> 

05 

• 

05 

o 

o 

CO 

X 

r- 

• 

CD 

o 

o 

CD 

05 

> 

00 

CD, 

o 

o 

* 

* 

z 

o 

h* 

< 

N 

N 

U 

L0 

CO 

X 

d 

c*. 

< 

a 

«- 

CO 

UJ 

o 

X 

Q. 

U_ 

polarization  difference  at  19.35  GHz.  Furthermore,  the  polarization 
difference  at  37  GHz  ia  less  sensitive  to  ice  type  than  either  the  37.0  GHz 
vertical  or  horizontal  brightness  temperatures.  Therefore,  the  difference 
Tgv(37)  -  TBjj(37)  can  be  employed  to  obtain  accurate  estimates  of  ice 
concentration. 

A  sensitivity  study  was  performed  to  determine  the  uncertainty  of  ice 
concentration  using  only  the  two  37  GHz  channels  (6,7).  For  winter 
conditions 

Tf  =  T  =  260  ±  10°K 

T  =  250  ±  10°K 
m 

T  =  271K 
w 

L  =  0.025  ±  0.025  cm, 

and  for  instrument  noise  contributions  of  AT  =  0  K  and  aT  =  1  K,  ice 
concentration  retrieval  accuracy  ranges  from  7%  when  C  =  100%  to  6%  for  low 
ice  concentrations.  The  uncertainty  for  low  ice  concentration  arises  mainly 
from  atmospheric  variability,  while  for  high  ice  concentration,  variability 
in  ice  emissivity  is  the  main  source  of  uncertainty.  These  analyses 
indicate  that  the  37  GHz  channels  will  provide  an  adequate  margin  for  ice 
concentration  measurements.  Thus,  the  proposed  algorithm  to  retrieve  ice 
concentration,  IC,  is 

IC  =  Aq  +  A1  CTbv(37)  -  Tbh(37)3,  (III. 48) 

where  Aq  and  A^  are  constants  computed  using  equations  III. 46  and  III. 47  for 
the  vertical  polarization  and  equivalent  equations  for  Tbb.  Climatologically 
significant  mean  values  of  ice  temperatures,  atmospheric  liquid  water 
content,  and  mean  values  of  ice  emissivities  are  used  as  input  to  these 
equations.  Aq  and  A^  vary  with  season.  They  are  presented  in  Section  IV. 

Dr.  Rene  Ramseier  of  the  Atmospheric  Environment  Service  of  Canada  has 
conducted  many  experiments  to  test  the  SSM/I  ice  algorithm  using  NIMBUS-7 
SMMR  37.0  GHz  data.  Slight  adjustments  of  coefficients  were  made  to 
accommodate  the  difference  in  satellite  geometry.  Results  from  this  study 
shows  good  agreement  in  ice  concentration  and  ice  edge  between  SSM/I  ice 


Ae  =  e  -  e.  =  0.05  ±  0.05 
v  h 

X  =  0.81  cm 
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’i 

p* 

■i 
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retrieval  and  ground  truth  data  collected  from  reconnaissance  flights,  and 
surface  observation. 

Ice  age  determination  (FY  or  MY)  can  be  derived  t-om  the  37  GHz 
vertical  brightness  temperature  value  and  computed  value  of  ice 
concentration  C.  This  is  accomplished  by  computing  the  effective  average 
ice  brightness  temperature  and  comparing  it  with  known  FY  and  MY  Tn  values. 

D 

Let  Tx  be  the  effective  average  ice  brightness  temperature  for  vertical 
polarization . 


e  1 

T  - -  [Tn„  -  T  ]  - —  <1-10  [f  T  +  (1  -  e  )T  ]. 

x  IC  BV  atm  I(,  w  w  w  s 


(III. 49) 


With  same  algebraic  manipulation,  this  reduces  to 


Tx  =  (Co  TBV  +  V/IC  "  C2*  (III. 53) 

If  Tx  is  greater  than  a  predetermined  value,  ice  age  is  designated  as  FY. 

The  coefficients  Cq,  C1  and  C 2  are  calculated  based  on  climatological 
values.  The  criteria  for  determination  of  ice  age  together  with  these 
coefficients  are  presented  in  Section  IV. 


IV.  RETRIEVAL  OF  GEOPHYSICAL  PARAMETERS 


Eleven  climates  are  defined  for  the  environmental  parameter  retrieval, 
for  both  the  land  and  the  ocean  cases.  These  climates  are  defined  in  Table 
IV. 1.  A  D  matrix  of  regression  coefficients  is  defined  for  each  of  the 
climates.  Note  that  there  are  four  transition  climates  within  the  eleven. 
The  D  matrices  corresponding  to  these  are  derived  from  the  neighboring 
climates.  For  the  sake  of  convenience,  code  numbers  are  designated  for  each 
of  these  climates  (see  Table  IV. 1). 

IV. A  Piecewise  Algorithm  Criteria 

The  piecewise  algorithm  employed  for  the  SSM/I  environmental  parameter 
retrieval  is  described  in  Subsection  II. C.  It  checks  certain  criteria  to 
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■ 

TABLE  IV. 1 

Definition  of  Climate  Codes 

• 

55 

Climate  Code  No. 

Definition 

Jc]  1 

Tropical-warm 

2 

Tropical-cool 

■l  3 

Lower  Latitude  Transition-warm 

M  4 

Lower  Latitude  Transition-cool 

^  5 

Mid.  Lat. -Spring/Fall 

6 

Mid.  Lat. -Summer 

4  7 

Mid.  Lat. -Winter 

y  8 

Upper  Lat.  Transition-cool 

9 

Upper  Lat.  Transition-cold 

ra  io 

Polar-cool 

*3  ii 

Polar-cold 

determine  which  parameters  to  retrieve  and  which  not  to.  It  is  designed  to 
avoid  meaningless  retrievals  and  to  use  computing  facilities  most 
effectively.  The  criteria  given  for  ocean  cases  include  the  critical 
temperature  of  the  19.35  H  brightness  temperature  ‘for  the  'maybe  rain'  case 
CMRO,  and  the  critical  brightness  temperature  difference  between  the  two 
polarizations  of  the  37  GHz  frequency  for  'maybe  rain',  CMRDO.  When  either 
of  the  two  criteria  is  met,  i.e.,  when  TD  ( 1 9H )  is  greater  than  CMRO  or  when 

D 

Tg(37V)  -  Tg(37H)  is  less  than  CMRDO,  the  algorithm  assumes  that  there  is 
rain.  The  criterion  for  'heavy  rain*  over  the  ocean,  CHRDO,  is  also  based 
on  the  difference  in  Tn  between  the  two  37.0  GHz  channels.  When  T  (37 V)  - 

D  D 

Tg(37H)  is  less  than  CHRDO,  the  'heavy  rain*  condition  is  assumed. 

For  land  retrievals,  there  is  a  criterion,  CFGL,  defined  to  test  for 
frozen  ground.  If  Tg(37V)  is  less  than  CFGL,  the  algorithm  assumes  frozen 
ground  and  considers  soil  moisture  retrieval  impossible.  The  retrieval 
algorithm  also  examines  the  37  GHz  brightness  temperatures  for  'heavy  rain' 
using  the  criteria  CHRL  and  CHRDL.  The  vertical  channel  brightness 
temperature  is  checked  against  the  criterion  CHRL.  If  it  is  less  than  CHRL 
and  the  difference  between  the  two  polarizations  is  less  than  CHRDL,  'heavy 
rain'  is  assumed  and  both  the  soil  moisture  and  the  cloud  liquid  water  over 
land  will  not  be  retrieved.  If  the  'heavy  rain'  condition  is  not  met,  the 
algorithm  then  checks  the  37  GHz  brightness  temperatures  using  criteria  CMRL 
and  CMRDL  in  a  similar  way  to  CHRL  and  CHRDL  for  the  'maybe  rain'  condition. 

If  no  rain  is  indicated,  the  rain  rate  over  land  and  the  liquid  water  over 

land  are  set  to  zero.  The  piecewise  algorithm  criteria  are  listed  in  Table  i 

IV. 2  for  the  eleven  climates. 

IV. B  Ocean  Retrievals 

IV.B.1  Ocean  Surface  Wind  Speed  (SW) 

The  response  of  brightness  temperature  at  19.35  H  &  V,  22.235  V,  and  37 
V  &  H  to  ocean  surface  wind  speed  is  demonstrated  in  Figure  IV. 1  for  a 
tropical  clear  atmosphere.  The  calculations  are  made  at  the  look  angle  of 
the  SSM/I.  Note  that  the  horizontal  channels  are  more  sensitive  to  wind 
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TABLE  IV. 2  Criteria  for  the  Piecewise  Algorithm 


Climate  Code  CMRO 


CMRO 

CMRDO 

CHRDO 

CFGL 

CHRL 

CHRDL 

CM 

190 

25 

10 

150 

9 

273 

10 

26 

190 

25 

10 

150 

273 

10 

26: 

190 

25 

10 

150 

273 

10 

26: 

190 

25 

10 

150 

273 

10 

26 

170 

25 

20 

240 

270 

10 

24i 

190 

25 

15 

150 

270 

10 

26 

160 

30 

20 

240 

270 

10 

241 

150 

35 

20 

240 

270 

10 

24i 

140 

35 

20 

270 

270 

— 

27i 

150 

35 

20 

240 

270 

10 

24i 

140 

35 

20 

270 

270 

— 

271 

CMRO 

CMRDO 

CHRDO 

CFO. 

CHRL 

CHRDL 

CHRL 

CHRDL 


Criterion  for  maybe  rain  over  ocean  for  T0 ( 1 9H ) 

Criterion  for  maybe  rain  over  ocean  for  TQ(37V)  -  Tn(37H) 

D  D 

Criterion  for  heavy  rain  over  ocean  for  Tg(37V)  -  Tg(37H) 
Criterion  for  frozen  ground  for  Tg(37V) 

Criterion  for  heavy  rain  over  land  for  Tg(37V) 

Criterion  for  heavy  rain  over  land  for  T0(37V)  -  T0(37H) 

Criterion  for  maybe  rain  over  land  for  Tg(37V) 

Criterion  for  maybe  rain  over  land  for  Tg(37V)  -  Tg(37H) 


for  heavy  rain  over  ocean  for  Tg(37V)  -  Tfi(37H) 
for  frozen  ground  for  Tg(37V) 


V£ 


ESS  TEMPERATURE  (DE 


speed  than  the  vertical  channels.  The  SSM/I  wind  speed  retrieval  algorithm 
is  given  by 


sw 

=  C  +  C  ,  T_ 

(19H)  +  C  , 

.  Td(22V)  * 

(IV. 1) 

sw,o  sw, 1  B 

sw. 

2  B 

+  Csw,3  TB 

(37 V)  +  C 

sw, 

4  Tb(37H), 

where  SM 

1  is  surface  wind  speed.  The  coefficients 

for  all 

eleven  climates 

are  listed  below: 

Climate 

Code  C 

C 

C 

C 

c  „ 

sw,o 

sw,  1 

sw,2 

sw,3 

sw,4 

1 

191.5600 

.4903 

-.4432 

-.9199 

.3577 

2 

168.3900 

.5366 

-.4548 

-.7656 

.2635 

3 

177.3150 

.3913 

-.2818 

-1.0083 

.4095 

4 

147.7600 

.5077 

-.3547 

-.7409 

.2333 

5 

127.1300 

.4788 

-.2546 

-.7162 

.2030 

6 

163.0700 

.2923 

-. 1204 

-1.0967 

.4612 

7 

95.9940 

.6106 

-.3034 

-.4638 

.0192 

8 

130.4200 

.3676 

-.1508 

-.8400 

.3056 

9 

117.5900 

.4225 

-.1899 

-.7096 

.2081 

10 

130.4200 

.3676 

-.1580 

-.8400 

.3056 

11 

117.5900 

.4225 

-.1899 

-.7096 

.2081 

IV. B. 2 

Precipitation  Over  Ocean  (RO) 

At 

lower  frequencies. 

such  as  19. 

35  GHz,  rain 

is  highly  absorptive  and 

results 

in  an  apparent  warm  brightness 

temperature 

over  a 

cold  background 

such  as 

the  ocean.  With  increasing  rain  rate,  the 

difference  in  brightness 

temperature  between  the  two  polarizations  at  a  given  frequency  tend  to 
decrease.  These  two  characteristics  are  demonstrated  in  Figure  IV. 2. 

The  SSM/I  data  channels  selected  for  retrieval  of  rain  over  ocean  are 
19.35H,  22.235V,  37.0  V  &  H.  The  algorithm  for  the  retrieval  is  given  by 


RO  =  C  +  C  T„ ( 1 9H )  +  C  -  T_(22V) 

ro  ,o  ro,  1  B  ro ,  d  b 


(IV. 2) 


*  Cro,3  TB(37V>  *  Cro,4TB<37H>- 


V-.  j 


35 


The  coefficents  for  all  climates  are 


Climate 

code 

Cro,o 

Cro,  1 

Cro,2 

.  Cro,3 

^ro,  4 

1 

210.2800 

.1217 

-.7829 

-.1830 

.0998 

2 

215.1800 

.1026 

-.8059 

-.1944 

.1354 

3 

173.0400 

.1938 

-.6500 

-.2291 

.0808 

4 

169.2900 

.1523 

-.6065 

-.3531 

.2162 

5 

123.4000 

.2019 

-.4070 

-.5117 

.2969 

6 

135.8000 

.2659 

-.5170 

-.2751 

.0618 

7 

114.5500 

.2708 

-.6228 

-.2836 

.2521 

8 

9.5432 

.1796 

-.2109 

.1214 

-.0753 

9 

24.1020 

.0825 

.1367 

-.3411 

.0843 

10 

9.5432 

.1796 

-.2109 

.1214 

-.0753 

11 

24.1020 

.0825 

.1367 

-.3411 

.0843 

IV. B. 3 

Cloud 

Water  (CWO) 

and  Liquid  Water  (LW0) 

Over  Ocean 

The  SSM/I  retrieval  algorithm  for  liquid  water  in  the  atmosphere 
includes  two  separate  components.  One  of  these,  CWO,  retrieves  the  amount 
of  liquid  water  contained  in  cloud  droplets,  defined  to  be  less  than  100  m 
in  diameter.  The  dominant  radiative  process  for  droplets  of  this  size  range 
is  absorption  for  the  SSM/I  frequencies.  The  response  at  19.35,  22.235  and 
37  GHz  to  integrated  cloud  water  over  a  mid-latitude  ocean  background  is 
shown  in  Figure  IV. 3.  The  data  channels  selected  for  the  retrieval  of  CWO 
are  the  same  as  for  all  the  other  ocean  parameters,  i.e.,  19.35  H,  22.235  V, 
37H  and  V. 


CWO  s  C  +  C 

CWO  ,o  CWO  ,1 


♦  C 


cwo,  3 


V19H)  +  Ccwo,2  V22V> 
V37V)  +  Ccwo,4  V37H)‘ 


(IV. 3) 
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BRIGHTNESS 


0 


Figure  IV. 3 


.05  .10  .15 

CLOUD  WATER  (gm/cm*) 

Brightness  temperature  vs.  cloud  water  over  mid-latitude 
ocean  at  19.  22  and  37  GHz 


The  coefficients  for  all  the  climate  codes  are  listed  below: 


Climate 

Code 

C 

CWO  ,0 

C 

CWO,  1 

C 

cwo, 2 

•  c 

cwo,  3 

r 

cwo, 4 

1 

-6.5240 

-.0159 

-.0044 

.0407 

.0045 

2 

-6.5292 

-.0151 

-.0046 

.0408 

.0041 

3 

-6.2070 

-.0141 

-.0038 

.0383 

.0035 

4 

-6.1372 

-.0094 

-.0054 

.0398 

-.0003 

5 

-5.7451 

-.0037 

-.0062 

.0387 

-.0047 

6 

-5.8903 

-.0123 

-.0031 

.0358 

.0025 

7 

-6.1411 

-.0007 

-.0061 

.0395 

-.0059 

8 

-4.0977 

-.0066 

-.0013 

.0230 

.0020 

9 

-3.6289 

-.0047 

-.0021 

.0212 

.0011 

10 

-4.0977 

-.0066 

-.0013 

.0230 

.0020 

11 

-3.6289 

-.0047 

-.0021 

.0212 

.0011 

The 

other  liquid  water 

retrieval 

algorithm 

for  the  ocean  is  desigi 

retrieve 

the 

integrated  liquid  water 

contained 

in  drops  of 

diameter  gr. 

than  100  m,  i.e.f  the  dimensions  of  rain  drops.  The  regression  equation 
for  the  retrieval  of  LWO  is 


0,0  =■  Cl-o,o  -  Cl».,  TB<,9H)  ♦  Clwo.2  V22V> 


(IV. 4) 


*  Clwo,3  V37V>  "  Clwo.4  V37H)- 


where  the  coefficient  for  the  eleven  climates  are  as  follows: 


Climate  Code 

c. 

lwo,o 

Clwo, 1 

Clwo,2 

Clwo,3 

Clwo, 

1 

50.7180 

.0300 

-.1881 

-.0493 

.0285 

2 

53.9140 

.0259 

-.2053 

-.4830 

.0374 

3 

40.7175 

.0405 

-.1353 

-.0768 

.0289 

4 

39.4580 

.0344 

-.1226 

-.0844 

.0342 

5 

25.0010 

.0429 

-.0399 

-.  1205 

.0310 

6 

30.7170 

.0510 

-.0824 

-.1043 

.0293 

7 

34.6250 

.0474 

-.0592 

-.1749 

.0623 

8 

.3091 

.0209 

-.0049 

-.0037* 

-.0142 

9 

-.8477 

.0048 

.0177 

-.0137 

-.0014 

10 

.3091 

.0209 

-.0049 

.0037* 

-.0142 

11 

-.8477 

.0048 

.0177 

-.0137 

-.0014 

*The  difference  in  sign  of  the  two  numbers  with  asterisks  is  apparently  a 
mistake.  The  Fleet  Numerical  Oceanography  Center  and  subsequently  Hughes 
Aircraft  Company  have  been  notified  of  this  error. 


IV. C  Land  Retrievals 

IV.C.1  Soil  Moisture  (SM) 

The  emissivity  of  soil  for  various  surface  soil  moisture  contents  is 
shown  in  Figure  IV. 4  for  the  1.43,  19  and  37  GHz  frequencies.  The  19.35  7  & 
H  GHz  channels  are  selected  for  the  retrieval  of  soil  moisture.  The  re¬ 
trieval  algorithm  for  soil  moisture  is  given  by 


SM  =  C  ♦  C  ,  T_(19V)  ♦  C  T_(19H). 
sm,0  3m, 1  B  sm,2  B 


The  coefficients  for  all  the  climate  codes  are  defined  as  follows: 


(IV. 5) 


Jf 

t 

w 

Climate  Code 

C 

sm  ,o 

Csm ,  1 

C 

sm. 

> 

1 

121.3000 

-0.4596 

0.0839 

2 

108.8100 

-0.3877 

0.0494 

,1, 

3 

109.5830 

-0.3872 

0.0414 

i 

4 

89.9820 

-0.2924 

0.0058 

». 

Y 

5 

71.1540 

-0.1970 

-0.0378 

l 

6 

97.8660 

-0.3148 

-0.0011 

i 

7 

85.8390 

-0.2818 

-0.0040 

8 

67.3730 

-0.1483 

-0.0874 

\ 

9 

0 

0 

0 

10 

63.5920 

-0.0995 

-0.1370 

3 

11 

0 

0 

0 

For  the  polar  cold  season  and  the  upper-latitude  cold  season  transition 
zones  (Climate  Code'  11  and  9),  land  is  regarded  as  frozen  and  soil  moisture 
is  not  retrieved. 

IV. C. 2  Precipitation  Over  Land  (RL) 

The  calculated  brightness  temperature  at  37  and  85.5  GHz  due  to  rain 
over  land  at  mid-latitudes  is  given  in  figure  IV. 5.  Rain  causes  a  reduction 
in  the  apparent  brightness  temperature  over  the  warm  land  background  because 
the  back  scatter  of  the  cold  upper  atmosphere  begins  to  dominate  the  forward 
scatter  from  the  land  surface  and  the  self  emission  from  the  atmosphere. 

This  effect  is  absent  over  the  ocean  due  to  the  much  lower  ocean  background 


40 


IVITY 


5  10  15  20  25 

PERCENT  MOISTURE 

Figure  IV. 4  Emissivity  at  1.43,  18  and  37  GHz  vs 
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Figure  IV. 5  Brightness  temperature  vs.  rain  rate  over  mid-latitude  land 
at  19,  37  and  85  GHz 


a  '  j  * 


and  the  coefficients  are  listed  below: 


Climate  Code 

Clwl,0 

Clwl,l 

Clwl  ,2 

1 

53.1520 

-0.1495 

• 

-0.0471 

2 

57.2120 

-0.1405 

-0.0713 

3 

51.3090 

-0.1580 

-0.0315 

4 

49.7870 

-0.0977 

-0.0873 

5 

42.3610 

-0.0548 

-0.1032 

6 

49.4650 

-0.1665 

-0.0158 

7 

52.2480 

-0.0687 

-0.1278 

8 

29.7660 

-0.0381 

-0.0725 

9 

0 

0 

0 

10 

17.1110 

-0.0213 

-0.0417 

11  0  0  0 
Since  liquid  water  does  not  generally  exist  in  polar  and  sub-polar  regions 
during  the  cold  season,  liquid  water  over  land  is  not  retrieved  for  climate 
codes  11  and  9* 

IV. C. 4  Cloud  Water  Over  Land  (CWL) 

Cloud  water  over  land  is  much  more  difficult  to  retrieve  than  cloud 
water  over  the  ocean  (6,7)  because  the  land  background  is  much  warmer  than 
the  ocean.  As  a  result,  at  the  lower  frequencies  there  is  little 
sensitivity  to  cloud  water.  Besides,  the  land  background  is  more  variable 
(7)  than  the  ocean.  The  brightness  variation  due  to  cloud  water  is  of  the 
same  order  of  magnitude  as  the  noise  level  of  the  land  surfaces.  On  the 
other  hand,  the  higher  frequency  channels  tend  to  be  so  sensitive  to  cloud 
water  that  they  serve  better  as  indicators  of  the  existence  of  clouds  than 
of  cloud  water  content  (6).  The  algorithm  developed  for  the  SSM/I  retrieval 
of  CWL  is  given  by 


The  coefficients  used  in 

retrieval 

for  all  the  climates 

are  defined 

below: 

Climate  Code 

C  . 
cwl  ,o 

Ccwl , 1 

C  ,  „ 
cwi,  a 

Ccwl,3 

Ccwl ,4 

1 

.6221 

.0036 

-.0031 

.0251 

-.0268 

2 

1.6545 

-.0008 

-.0006 

.0362 

-.0384 

3 

.8821 

.0006 

-.0007 

.0271 

-.0283 

4 

1.5350 

.0014 

-.0030 

.0362 

-.0377 

5 

1.4154 

.0036 

-.0054 

.0362 

-.0369 

6 

1.1421 

-.0024 

.0017 

.0290 

-.0298 

7 

1.2078 

-.0050 

.0033 

.0415 

-.0401 

8 

1.5366 

.0002 

-.0019 

.0328 

-.0342 

9 

.9329 

-.0045 

.0028 

.0351 

-.0333 

10 

1.6578 

-.0033 

.0017 

.0294 

-.0314 

11 

.6579 

-.0040 

.0023 

.0286 

-.0265 

IV. D  Ice  Retrievals 


The  SSM/I  retrieval  algorithm  for  ice  is  derived  differently  than  those 
for  ocean  and  land  retrievals.  Certain  relationship  among  the  brightness 
temperatures  and  ice  parameters  were  predetermined  (See  Equations 
III. 48-50),  and  coefficients  of  these  relationship  are  derived  through 
analytic  expressions  and  climatological  data  as  described  in  subsection 
III.C.3.  The  coefficients  AQ  and  A1  from  equation  III. 4. 8  and  Cq,  Cn 
from  equation  III. 50  for  the  derivation  of  ice  concentration,  IC,  and  the 
effective  average  ice  brightness  temperature,  Tx,  are  listed  below: 


Polar-cool 

Polar-cold 


‘1 


'1 


T  (MY) 
c 


1.1880 

1.1780 


-0.0176 

-0.0176 


-219.0600 

-217.7000 


1.1050 

1 . 1050 


193.6400 

192.2800 


228.00 

218.00 


where  T  (MY)  is  the  critical  temperature  for  the  determination  of  ice  type, 
c 

If  T„  is  greater  than  Tc(mY),  the  ice  is  defined  to  be  first  year,  otherwise 
it  is  designated  to  be  MY  ice. 

Cloud  water  over  ice  (CWI)  is  retrieved  by  the  SSM/I  algorithm  even 
though  it  is  not  required  by  the  specifications.  The  data  channels  employed 
are  19.35  V  A  H  and  37.0  V  4  H. 


3 


o 


A 
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(IV. 9) 


CWI  =  c„ ^  „  +  C  .  .  T_(19V)  +  C  .  . 
cwl  ,o  cwi,1  B  cwit2 


*  Ccwl,3  VOT>  *  W.4 


The  coefficients  employed  are  listed  below. 


TgC 19H) 
Tb(37H). 


cwi  ,o 

cwi,  1 

cwi  ,2 

cwi,  3 

cwi ,  4 

Polar-cool 

.3415 

-.0010 

-.0048 

-.0667 

.0744 

Polar-cold 

-.2082 

-.0003 

-.0121 

-.0320 

.0468. 
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